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ABSTRACT 


The  operation  of  a  strapped-down  inertial  navigation 
system  in  a  high  dynamic  environment  was  simulated.  Its 
performance  was  assessed  using  a  laser  gyro  and  a  quartz- 
flexure  accelerometer  as  the  sensors.  Using  deterministic 
values  for  the  sensor  errors,  the  effects  of  the  sensor 
errors  on  the  INS  performance  was  assessed.  Emphasis  is; 
placed  on  the  effects  that  the  sensor  errors  have  on  the  INS 
errors  in  a  high  dynamic  environment.  Structural  mode 
effects  on  INS  performance  are  also  evaluated.  Following 
this  an -analysis  of  the  sensitivity  of  the  INS  errors  in 
position  and  velocity  to  each  inertial  sensor  is  determined. 
Using  this  sensitivity  analysis,  the  sensor  paramenters  re¬ 
quired  to  obtain  a  performance  specification  of  -  1  ni^/hr 
position  error  and  ^'3  ft/sec  velocity  error  are  defined. 
Finally,  recommendations  of  areas  for  improvement  are  present¬ 
ed. 
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CHAPTER  I 


INTRODUCTION 


Background 

It  would  be  difficult,  if  not  impossible,  to  determine 
when  the  first  inertial  navigation  system  was  developed. 

The  first  patent  for  an  assemblage  of  three  gyros  and  a 
gravity  vertical  as  a  navigation  system  was  given  in  1924 
to  an  American,  C.G.  Abbot.  Since  that  time  much  effort 
has  been  spent  in  regards  to  inertial  navigation  systems 
(INS).  During  the  1930's  the  concept  of  using  accelero¬ 
meters  mounted  on  a  stabilized  platform  was  utilized  (Ref 
5«8-12). 

Through  the  years  many  different  types  of  INS  have  been 
designed,  but  they  can  be  classified  under  two  distinct 
groups,  strapped-down  and  gimbaled  systems.  The  distin¬ 
guishing  feature  of  these  two  groups  is  the  way  in  which 
their  orientation  with  respect  to  inertial  space  is  main¬ 
tained.  In  the  case  of  most  gimbaled  systems  the  platform 
is  mounted  with  gimbales  such  that  it  remains  fixed  in  some 
orientation  with  respect  to  inertial  space  regardless  of 
the  movements  of  the  body  on  which  it  is  mounted.  The 
strapped-down  systems  are  fixed  to  the  body  and,  therefore, 
the  relation  of  the  platform  frame  to  inertial  space  is 
maintained  by  updating  the  transformation  matrices  using 
the  information  received  from  the  gyros.  Because  of  the 


1 


•xtra  computational  load  required  for  the  strapped-down  sys- 
to  maintain  orientation  with  inertial  space,  this  group  of 
systems  received  little  attention  until  the  1960's.  As  com¬ 
puter  technology  progressed,  the  size  and  weight  of  equip¬ 
ment  necessary  to  handle  this  additional  computational  load 
was  drastically  reduced  and  more  emphasis  was  placed  on 
strapped-down  systems.  As  further  refinements  have  been 
and  are  being  made  in  the  area  of  micro  computers  during 
the  seventies,  development  of  the  laser  gyro  also  pro¬ 
gressed.  With  these  advances  in  technology,  the  potential 
of  the  strapped-down  systems  have  grown.  Possible  advan¬ 
tages  that  strapped-down  systems  could  have  over  gimbaled 
systems  include  smaller  size,  less  weight,  and  greater 
maintainability  and  reliability.  With  these  factors  in 
mind,  future  commitments  of  resources  and  money  for  fur¬ 
thering  the  development  of  the  strapped-down  INS  is  desira¬ 
ble  (Ref  2 1 1 ) . 

Statement  of  Problem 

Past  analysis  of  inertial  navigation  systems  has  shown 
that  the  gimbaled  systems  are  more  accurate  than  the  strapp- 
down  systems.  The  implications  are  that  the  reason  for  this 
difference  in  performance  of  the  strapped-down  and  gimbaled 
INS  is  due  to  the  harsher  environment  that  the  strapped-down 
is  subjected  to.  An  analysis  of  the  error  propogation  of 
the  sensor  errors  should  give  insight  into  the  effect  of  the 
sensor  errors  on  the  total  system  errors  (Ref  7 » 37-38).  This 
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would  then  be  a  basis  for  working  at  improving  the  strapp- 
'  down  INS  performance.  In  undertaking  this  error  analysis, 

a  digital  strapped-down  INS  simulation  recently  developed 
by  the  Charles  Stark  Draper  Lab  under  contract  to  the  Air 
Force  Avionics  Lab  at  Wright-Patterson  AFB  is  utilised. 

Plan  of  Attack 

Using  the  strapped-down  INS  simulation  developed  by  the 
Draper  Lab,  a  series  of  INS  error  time  histories  will  be 
generated.  From  examining  the  time  histories,  the  effect 
of  the  sensor  errors  on  the  INS  errors  will  be  evaluated. 

By  changing  the  values  of  the  sensors  errors,  the  sensiti¬ 
vity  of  INS  errors  to  sensor  errors  will  be  developed. 
Finally,  the  sensor  parameters  necessary  to  achieve  a  posi- 
^  tion  error  of  «  1  nm/hr  and  a  velocity  error  of  ^3  ft/sec 

will  be  determined. 

Sequence  of  Events 

A  general  overview  of  the  capabilities  of  the  strapped- 
down  INS  simulation  is  presented  first.  Following  this,  a 
detailed  description  of  the  sensor  and  navigation  models  is 
presented.  Also,  a  description  of  the  flight  profile  and 
sensor  parameters  used  in  the  simulations  is  given.  Based 
on  the  simulation  analysis  of  the  effect  of  sensor  errors 
on  INS  errors  is  performed.  Conclusions  are  then  drawn 
from  this  analysis. 

(. 
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CHAPTER  II 


Description  of  Computer  Simulation  and  Navigation  System 
and  Sensor  Model  and  Their  Implementation  in  Obtaining 
Data  for  an  Error  Analysis 

The  strapped-down  inertial  navigation  system  simulation 
referenced  in  the  introduction  provided  the  data  that  is 
used  in  the  error  analysis  that  is  presented  in  later 
chapters.  The  simulation  is  a  digital  computer  program 
coded  in  FORTRAN  and  is  implemented  on  a  GDC  6600  computer. 
The  program  was  designed  for  ease  of  user  interaction,  and 
proved  to  be  a  useful  tool  for  analyzing  inertial  navigation 
system  errors. 

This  section  provides  a  brief  introduction  to  the  capa¬ 
bilities  of  the  strapped-down  simulation.  A  detailed  des¬ 
cription  of  the  navigation  system  and  its  associated  iner¬ 
tial  sensors  are  then  presented.  Finally,  the  flight  pro¬ 
file  used  to  generate  the  simulations  for  the  error  analysis 
is  deadribed. 

Overview  of  Strapped-Down  Simualtion  Capabilities 

The  strapped-down  inertial  navigation  system  simulation 
(SDINS)  employs  an  F-4  aircraft  model.  Nonlinear  six-de- 
ree-of-freedom  equations  of  motion  and  nonlinear  aero  dyna¬ 
mics  valid  beyond  the  stall-angle-of -attack  are  implemented. 
Bending  effects  as  well  as  three  longitudinal-and  three 
lateral -directional  structural  modes  are  included. 
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The  capability  for  perturbation  inputs  such  as  winds  and 
gusts  is  also  available.  The  other  models  included  are 
those  for  a  flight -control  system,  navigation  system,  and 
their  associated  sensors,  as  well  as  a  gravity  model.  In 
the  work  done  for  this  paper,  a  rotating  spherical  earth 
model  is  used  for  the  computations,  but  an  elliptical  model 
is  provided  in  the  simulation. 

This  is  only  a  brief  overview  of  the  strapped-down  INS 
simulation  program.  Since  this  work  is  only  concerned  with 
the  details  of  the  actual  navigation  model  and  its  associat¬ 
ed  sensors,  this  protion  of  the  system  is  now  described  in 
detail. 

Navigation  System  Model 

Presented  here  is  the  algorithm  used  in  the  strapped- 
down  INS  simulation  to  represent  a  strapped-down  Local 
Vertical  Wander  Azimuth  (LVWA)  navigation  system.  The 
following  discussion  of  this  algorithm  is  based  on  infor¬ 
mation  taken  form  reference  (2),  The  navigation  model  can 
be  divided  into  two  major  subseotions,  a  low-speed  loop 
computed  at  5  hertz  and  a  high-speed  loop  computed  at  50 
hertz.  This  representation  results  form  the  need  to  up¬ 
date  the  body-to-computational  transformation  matrix  at  a 
high  frequency  in  order  to  maintain  accuracy »  while  the 
other  computations  need  not  be  performed  at  this  high  of 
an  iteration  rate. 


The  high  speed  navigation  loop  consists  of  updating  the 
body-to-computational  or  local -vertical  wander  azimuth 
transformation  matrix ,  and  transforming  the  navigation  sen¬ 
sor  velocities  from  the  body  frame  to  the  computational 
frame  and  the  summation  of  these  velocities  over  the  low 
speed  computational  time  period.  The  generation  of  the  bo¬ 
dy-to-computational  transformation  matrix  and  transforma¬ 
tion  of  the  velocities  are  staggered  in  timei  with  the 
transformation  matrix  being  computed  midway  between  the 
times  at  which  the  transformation  of  the  velocities  occur. 
The  update  of  this  transformation  matrix  is  achieved 
using  a  third  order  quaternion.  It  is  from  this  body- 
to-computational  frame  transformation  matrix  that  the 
attitude  angles  are  computed  for  comparison  with  the  cor¬ 
responding  true  angles  generated  in  the  vehicle  portion  of 
the  simulation. 

The  function  of  the  low-speed  portion  of  the  system 
model  is  to  perform  the  necessary  operations  on  the  sen¬ 
sor  signals  in  order  to  find  position  and  velocities.  All 
signals  are  defined  in  the  computational  frame,  but  a  trans 
formation  matrix  for  transforming  from  computational  frame 
to  earth  frame  is  maintained.  It  is  from  this  transforma¬ 
tion  matrix  that  latitude,  longitude,  and  wander  angle  are 
obtained  (Appendix  C  derives  the  transformation  in  terms  of 
these  quantities).  The  velocities  are  taken  from  the  inte¬ 
gration  of  the  accelerometer  signals.  A  further  integra¬ 
tion  of  the  vetical  channel  velocity  yields  altitude. 
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Since  the  vertical  channel  of  an  inertial  navigation  system 
is  unstable  without  aiding,  a  third-order  damping  scheme  is 
applied. 

The  equations  implemented  in  the  navigation  model  are 
now  examined.  Appendix  A  provides  a  review  of  the  notation 
and  reference  frames  that  are  used.  It  should  be  noted 
that  the  local  vertical  wander  azimuth  frame  is  the  same 
frame  as  the  computational  frame  and  it  is  referred  to  as 
the  computational  frame  in  the  following  discussion.  The 
x,  y,  and  z  coordinates  of  this  frame  are  defined  in  an  up, 
east,  and  north  sense  (Ref  Ii49-51). 

Navigation  Model  Equations 

Figure  1  shows  a  block  diagram  of  the  basic  inertial 
navigation  system  model.  The  accelerometers  provide  a 
measure  of  the  specific  force.  The  measurement  is  in  the 
body  frame  and  must  be  transformed  to  the  computational 
frame .  Thus , 


*  c£  £b  (1) 

where  £  represents  specific  force,  and  c£  is  the  trans¬ 
formation  matrix  to  transform  from  the  body  frame  to  the 
computational  frame,  and  c  and  b  represent  computational 
frame  and  body  frame,  respectively. 

The  block  labeled  I  represents  the  integration  of  the 
velocity  differential  equation  that  is  derived  in  Appendix 
B.  This  equation  is 


? 
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Gyros 


ic  ■  lc  ♦  £c  -  (2 w°e  ♦  ^c)  X  vc 


(2) 


where  g  represents  the  acceleration  due  to  gravity,  wie 
represents  the  rotation  rate  of  the  earth  frame  with 
respect  to  the  inertial  frame,  and  w^c  represents  the 
rotation  rate  of  the  computational  frame  with  respect  to 
the  earth  frame. 

By  dividing  the  proper  components  of  V?  the  compo¬ 
nents  of  wc  are  computed.  (The  block  labeled  w  repre- 
0  c 

sents  these  calculations).  They  are 


-w‘v  »> 


(3) 


and 


<ct  ■  vec/(Ro  *  h>  <4> 

where  R0  is  the  radius  of  the  earth  (assuming  spherical 

earth)  and  h  is  the  altitude  above  the  earth's  surface. 

(w  «0  since  the  x  axis  is  the  vertical  axis).  is 

ecx  ^oc 

then  fed  back  to  be  used  in  the  velocity  differential 
equations.  is  also  used  to  forraualte  C®,  the  trans¬ 

formation  matrix  to  transform  a  vector  from  the  earth 
frame  to  the  computational  frame.  To  do  this,  (using  the 

fact  wc  ■  -wc  ),  the  skew  symmetric  form  of  -w?a  is  found, 
ec  ce  v  e 
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-w_ 


C6t> 


wcet 

wcey 


cex 


-w 


cex 


(5). 


Then  this  matrix  is  used  in  the  formulation  of  the  rate 
of  change  of  a  direction  cosine  matrix i 


wck  rc 
ce  e 


(6) 


It  is  then  from  Cc  that  latitude,  longitude,  and  wander 
angle  are  calculated.  (Appendix  C  provides  the  deriva¬ 
tion  of  in  terms  of  these  quantities). 

The  gyros  provide  a  measure  of  the  angular  rate  of 
the  body  frame  with  respect  to  inertial  space,  w^.  This 
is  then  differenced  with  w^c  resulting  in 


b  b 
-ib  "  -ic 


(7) 


where  w^,.  the  rotation  rate  of  the  computational  frame 
relative  to  inertial  space  is  determined  by  adding  w£e 
"to  and  then  transforming  to  the  body  frame. 


(8) 


The  quantity  the  angular  velocity  of  the  earth  frame 

relative  to  inertial  space,  is  used  in  the  velocity  differ¬ 
ential  equation  as  well  as  in  equation  (8).  Assuming  a 
spherical  earth  rotating  at  a  constant  rate, 
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0 


0 


(9) 


where  w^e  is  the  rate  of  the  earth's  rotation,  15°/hr. 
Thus , 


Cc  w? 

e  —  le 


(10) 


Completing  the  gyro  loop  of  the  diagram  in  figure  1  re¬ 
quires  the  calculation  of  c£.  In  order  to  maintain  the 
desired  accuracy  in  calcualting  c£,  a  third  order  quater¬ 
nion  update  is  used  instead  of  the  direction  cosine  dif¬ 
ferential  equation.  However,  just  as  wf  is  needed  to 
c  b 

calculate  Ce,  wfcc  is  needed  to  perform  the  quaternion 
update.  Following  the  calculation  of  c£,  the  Euler  angles, 
roll,  pitch,  and  yaw,  are  computed  using  the  elements  of 
c£  (Ref  2i22-3l).  (Appendix  C  provides  a  derivation  of 
c£  in  terms  of  the  Euler  angles). 

6x=-g0  [l-2(h/R0)+3(h/RQ)2]  ♦  w?e*R0(l+h/R0)  (l-sin2L)  ( 11 ) 
gy=-w2e,Ro(l-«‘h/R0)  sinL  (-cosL  sin  d  )  (12) 

gz= ”wie * Ro ( 1 +h/^o ^  sinL  (°osL  cos  *  )  (13) 

o 

where  gQ  is  the  gravity  constant  (gQ=32.1725  ft/sec  ), 

Rq  is  the  radius  of  the  earth  (RQ=20,860,000  feet),  L 
is  the  latitude,  and  JL  is  the  wander  angle. 
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X 


Figure  2.  Orientation  of  Aircraft  Body  Axis. 

A  final  note  on  the  navigation  algorithm  before  examin¬ 
ing  the  sensor  models  that  are  used  to  generate  the  naviga¬ 
tion  algorithm  inputs  concerns  the  generation  of  errors  in 
position,  velocity,  and  attitude.  These  errors  are  generat¬ 
ed  by  the  differencing  the  navigation  model  outputs  with 
those  from  the  aircraft  equations-of-motion.  Therefore,  the 
error  quantities  include  effects  from  computational  errors 
(Ref  1»  58). 


Sensor  Model 

Figure  2  shows  the  orientation  of  the  aircraft  body 
axis.  There  are  three  gyros  and  three  accelerometers  model¬ 
ed  for  the  navigation  system.  Each  gyro  measures  the  ro¬ 
tation  rate  about  one  of  the  principle  body  axis.  Similarly, 
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there  is  an  accelerometer  located  on  each  body  axis.  It  is 
with  reference  to  this  configuration  that  the  discussion  on 
the  gyro  and  accelerometer  models  is  based.  Also,  only  the 
roll-rate  (Z  axis)  gyro  model  and  the  lateral  (Y  axis)  ac¬ 
celerometer  model  are  presented  since  they  are  representative 
of  the  others.  The  accelerometer  model  is  presented  first 
and  is  followed  by  the  gyro  model. 

Accelerometer  Model 

The  accelerometer  implemented  in  the  strapped-down  INS 
simulation  is  representative  of  the  quartz-flexure  class  of 
accelerometers.  The  accelerometer  error  sources  modeled  in¬ 
clude  misalignment,  bias,  scale  factor,  cross-coupling,  and 
nonlinearity.  The  model  assumes  that  the  statistical  data 
on  these  error  sources  for  a  specific  instrument  are  sup¬ 
plied  as  a  simulation  input.  Also,  the  model  has  a  limiter 
to  limit  the  accelerations  to  *  10 >g. 

Figure  3  shows  a  block  diagram  of  the  lateral  acceler¬ 
ometer  model.  The  input  signal,  A,  represents  the  rigid 
body  acceleration  including  bending  effects.  This  signal  is 
integrated  using  a  rectangular  integration  scheme.  Thus, 

AyI  ■  AyR  4  V*-11  .  it  (i4> 

2 

where  the  subscript  I  represents  an  integrated  signal, k  and 
k-1  represents  the  iteration  number,  and  At  is  the  period  be¬ 
tween  samples.  At  this  point,  the  effects  of  the  misalignment 
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error  are  accounted  for.  The  coefficient  of  misalignment, 
y  ,  is  multiplied  by  the  respective  signal  and  then  added 
to  the  integrated  lateral  acceleration.  Therefore, 

AyIM  “  Ayl  *  Axl  +  A  Azl  (15) 

where  the  subscript  M  represents  signal  with  misalignment 
effects  and  represents  the  other  accelration  signals  after 
they  are  integrated.  This  signal  is  now  fed  through  a  first 
order  lag  to  represent  the  delays  of  the  instrument. 

AyIML  '  AyIM<l/Wstl))  <l6) 

where  O'  =  time  constant  of  the  lag  and  subscript  L  indicates 
the  signal  passes  through  a  lag.  The  bias  error  is  added  to 
the  accelerometer  signal  after  it  passes  through  the  lag  net¬ 
work. 


AyIMLB  =  Ay XML  *  6  BIAS 


At 


(17) 


where  subscript  B  represents  signal  with  bias  effects  and 
€  bias  rePresents  the  statistical  value  for  the  bias  error 
as  stated  in  the  specifications  for  the  accelerometer  used. 
Before  the  quantization,  the  nonlinearity  error  and  cross¬ 
coupling  error  are  accounted  for, 


AyIMLBNC*  AyIMLB  *  (CyxAxIMLB  *  CyyAyIMLB)AyIMLB 


(18) 
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Table  Z 


Quartz-I-'Icxure  Accelerometer  Nominal  Values  (Ref  1»70) 


Parameter 

Symbol 

Value 

Misalignment  Coefficients 

Time  Constant  of  Lag 

Scale  Factor 

Bias  Error 

Cross-Coupling  Coefficients 

Nonlinearity  Coefficients 

Scale  Factor  Error 

h'h 

r 

sfa 

4  . 

Bias 

cyx 

cyy 

5  x  10"5  rad 

1  x  10'^  sec 
32,000  pps/g 

50  g 

19-6  x  10"5  rad/sec 
40  g/g2 

100  ppm 

where  subscripts  M  and  C  represent  nonlinearity  and  cross¬ 
coupling  coefficient  and  C  is  the  nonlinearity  coefficient. 
This  signal  is  now  converted  to  pulses  by  the  scale  factor, 
SFa,  and  then  quantized. 

Table  1  shows  a  listing  of  the  parameter  values  that 
are  used  for  the  accelerometer  in  the  simulations  that  were 
made  for  this  work.  The  error  parameter  values  are  the  de¬ 
terministic  tolerances  as  listed  in  manufacturer’s  perfor¬ 
mance  specifications  for  a  quartz-flexure  accelerometer. 

These  figures  are  representative  of  presently  available 
instruments,  but  do  not  necessarily  reflect  the  limit  of 
achievable  performance. 
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The  quantization  takes  an  integer  count  of  the  value 
of  the  present  signal  plus  the  value  of  the  fractional  part 
of  the  previous  sample  period  integer  count.  This  integer 
count  is  then  divided  by  the  scale  factor  and  the  effects 
of  the  scale  factor  error  are  accounted  for  resulting  in 
the  incremental  value  of  the  lateral  velocity, 

*y<l/SPA)(1’EA)AyIMMNM  (19) 

where  EA  is  the  scale  factor  error  and  subscript  Q  repre¬ 
sents  the  signal  including  quantization  and  scale  factor. 
The  incremental  velocities  are  then  summed  up  in  the  na¬ 
vigation  system  subroutine  of  the  strapped -down  INS  sim¬ 
ulation  to  be  used  in  the  calculation  of  velocity.  This 
summation,  as  was  mentioned  in  the  navigation  algorithm 
description,  is  done  at  a  higher  iteration  rate  than  the 
rate  at  which  the  velocity  differential  equation  is  solved 
(50  hertz  as  opposed  to  5  hertz) (Ref  1i47-49). 

Gvro  Model 

Two  types  of  gyros  are  modeled  in  the  strapped-down 
INS  simulation.  They  are  laser  gyro  and  two-degree-of- 
freedom  mechanical  gyro.  There  are  three  of  each  type 
modeled,  one  for  measuring  rotation  rate  about  each  of  the 
three  body  axis.  The  roll-rate  laser  gyro  is  described 
here  and  then  differences  of  the  two-degree-of-freedom 
gyro  are  pointed  out.  Once  again,  statistical  data  on 
the  instrument  error  is  assumed  available.  Error  sources 
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modeled  include  bias,  misalignment,  and  scale  factor. 

An  illustration  of  the  block  diagram  for  the  roll  rate 

laser  gyro  is  shown  in  Figure  4.  w„  the  angular  inertial 

«  » 

velocity  about  the  body  roll  axis  including  bending  effects 
is  fed  into  an  intrator.  Once  again  a  rectangular  integra¬ 
tion  routine  is  used.  Therefore, 


£ 


zl 


wzk  *  wz(k-l ) 


At 


(20) 


After  interpretation,  the  effects  of  misalignment  error 
are  accounted  for,  resulting  in 


•^zIM  cAl  +  A  Acl  (21) 

where and  /ii  represent  the  coefficients  of  misalign¬ 
ment  for  their  respective  signals.  The  bias  error  is  now 
added  to  the  roll  rate  signal.  It  is  a  function  of  the 
time  period  and  not  of  a  rotation  magnitude.  Also,  at 
this  point,  the  random  walk  is  accounted  for,  resulting 
in, 


^ zim  "  Aim  +  €bias  *  RL«(  at)*  (22) 

where  is  the  bias  error  and  R^  is  the  random  walk 

error.  ^RXi"ErtALK*  where  EWALK  is  'tolerance  spe- 

cifioatiDn  and  RAN  is  a  random  number).  The  signal  is 

now  converted  to  pulses  via  the  scale  factor,  SF^i 


Ring  Laser  Roll  Rate  Gyro  Model  (Ref  Ii46) 


^zIMBS  *  ^bIMB^SFG^  *23* 

This  signal  is  then  quantized  as  in  the  case  of  the  accel- 
eroroter  model,  and  then  is  multiplied  by  the  inverse  of  the 
scale  factor.  The  error  in  scale  factor  is  then  taken  into 
account  resulting  in  incremental  value  of  the  roll  axis 
rotation. 


A©a  a  ^1/SFG^AlMB3^1+ 

where  is  the  error  in  scale  factor.  This  signal  is  now 
summed  in  the  high  speed  loop  of  the  navigation  system  mo¬ 
dels  before  being  used  in  the  calculation  of  the  velocity 
and  position  (Ref  lt46-4?). 

Figure  5  shows  a  block  diagram  of  the  roll  channel  of 
the  two-degree-of-freedom  mechanical  gyro  model.  This 
model  is  the  same  as  the  laser  gyro  except  that  it  ac- 
counts  for  errors  peculiar  to  the  two-degree-of-freedom 
mechanical  gyro.  These  errors  are  the  acceleration  depen¬ 
dent  and  angular  dependent  terms.  Also,  the  laser  gyro  model 
requires  a  lag  to  account  for  delays  that  are  not  present 
in  the  laser  gyro  model  (Ref  3 i 22-25). 

The  values  of  the  parameters  used  in  the  laser  gyro 
and  two-degree-of-freedom  mechanical  gyro,  when  obtaining 
the  simulations  for  the  error  analysis,  are  listed  in  Tables 
II  and  III.  These  values  are  representative  of  presently 

i 

available  sensors. 


20 


Table  II 


Laser  Gyro  Nominal  Values  (Kef  1»?0) 


Parameter 

Symbol 

Value 

Wisal ignmont  Coefficients 

■4, 

5  x  10”^  rad 

Scale  Factor 

sfg 

1.57  sec/ pulse 

Bias  Fri  er 

Bias 

0.01  °/hr 

Scale  Factor  Hrror 

£-g 

5  x  10” ^  ppm 

Fandom  Walk  iirror 

^  walk 

0.005  °/hr 

Since  the  two  types  of  gyro  models  are  available  for 
the  simualtions,  a  comparison  of  the  results  of  the  simu¬ 
lation  ..using  each  gyro  model  is  made.  Figure  6  shows  the 
INS  error  time  histories  obtained  for  both  simulations. 

The  laser  gyro  provides  better  overall  performance  than  the 
two-degree-of-freedora  mechanical  gyroj  and  it  has  the  addi¬ 
tional  errors  due  to  acceleration  and  angular  rate  dependent 
terms.  Based  on  these  results »  the  laser  gyro  is  chosen 
for  the  simulations  used  in  this  work. 
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TABLE  III 


Two-Degree-of-Freedom  Gyro  Nominal  Values  (Ref  306) 


Parameters 

Symbol 

Misalignment  Coefficients 

Scale  Factor 

SPG 

Bias  Error 

Scale  Factor  Error 

^•Bias 
£  TDK 

g-Dependent  Errors 

fin 

£  -Dependent  Errors 

TDF4 

TDF5 

tdf6 

g^-Dependent  Errors 

TDF? 

tdf8 

tdf9 

w  -Dependent  Errors 

aW 

TDFi0 

TDF11 

TDF12 

Z-Output 

-Axis 


10'4  rad 

1*57  see 
pulse 
0.01  dcp. 

'hr 

50  ppm 


0.04  £)££ 

0.0!  I 


0.00  sU:fi 
0.005  £2 


0 . 04  fc.p 

c 

0,00  deg 
80  h  2 

20  (ra» 

see 


Y-Output 

-Axis 


0.01 

0.04  deg 

0.04  &2 
a 

0.00 

20 

80 


of  Position  Error  Tine  Histories  with  only  Gyro  Errors  Present 
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As  is  eluded  ta  in  the  introduction,  part  of  the  errors 
in  a  strupped-down  INS  result  form  the  high  dynumic  envi¬ 
ronment  that  the  sensors  are  subjected  to  when  the  navi¬ 
gation  system  is  used  on  a  fighter-type  aircraft.  Time 
histories  for  the  INS  errors  in  position  velocity  and  at¬ 
titude  are  generated  using  the  strapped -down  INS  simula¬ 
tions  so  that  the  effects  of  the  high  dynamic  environment  - 
can  be  assessed.  Therefore,  a  flight  profile  containing 
several  maneuvers  typical  of  a  tactical  aircraft  mission 
is  chosen.  Figure  7  shows  the  flight  profile  implemented 
for  this  work.  The  initial  start  of  the  simulation  assumes 
that  the  aircraft  is  in  straight  and  level  flight  at  a 
cruise  speed  of  mach  .5  at  an  altitude  of  5000  feet,  and 
a  heading  of  north.  At  100  seconds  into  the  simulation, 
the  aircraft  descends  to  500  feet  while  increasing  its 
velocity  to  mach  .75.  After  completion  of  the  descent, 
the  aircraft  performs  a  snap  roll.  Following  the  roll,  at 
about  360  seconds  into  the  flight,  a  weapon  delivery  is  sim¬ 
ulated.  This  maneuver  consists  of  a  dive  at  a  12.5  degree 
angle  followed  by  a  3*5g  pull-up.  After  a  leveling  off 
from  the  weapon  delivery,  the  aircraft  performs  another 
snap  roll  followed  by  an  evasive  maneuver.  The  evasive  ma-' 
euver  consists  of  a  climb  at  a  flight  path  angle  of  20°  with 
a  bank  of  45°  to  the  right  for  seven  seconds  foiblowod  by  a 
bank  to  the  loft  of  45°  for  seven  seconds,  leveling  off  at 
an  altitude  of  5000  feot  and  a  speed  of  mach  .75.  Upon 
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Figure  7«  Description  of  Flight  Profile  (Ref  4iH73) 
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leveling  off  at  5000  feet,  the  aircraft  performs  a  third 
snap  roll  followed  by  a  climb  to  14,500  feet  during  which 
its  speed  is  decreased  to  mach  .5*  The  final  maneuver  sim¬ 
ulated  is  a  loiter.  During  this  phase  of  the  flight,  the 
aircraft  turns  through  three  revolutions  at  a  30°  bank 
before  returning  to  straight  and  level  flight  on  a  north 
heading.  This  maneuver  is  completed  at  about  the  1620  se¬ 
cond  time  point.  The  aircraft  then  maintains  straight  and 
level  flight  for  the  remainder  of  the  flight  time  (Ref  4i 
133-135). 

Summary 

This  chapter  has  provided  a  description  of  the  models 
implemented  and  assumptions  made  in  using  the  strapped- 
down  simulation  to  obtain  the  INS  error  time  histories 
used  in  the  deterministic  error  analysis  that  is  to  fol¬ 
low.  This  is  not  meant  as  a  complete  description  of  the 
strapped -down  simulation  and  does  not  reflect  the  full 
capabilities  of  this  computer  simulation. 


CHAPTER  III 


Effects  of  the  Sensor  Errors 
and  Struc t ural  Modes  on  the  INS  Errors 


This  chapter  presents  a  deterministic  evaluation  of  the 
effects  the  sensor  errors  have  on  the  errors  of  the  strapped 
down  INS.  Time  histories  of  the  INS  errors  in  attitude • 
velocity i  and  position  are  obtained  using  the  strapped- 
down  INS  simulation.  The  simulations  are  made  for  the  case 
with  each  sensor  error  source  isolated  as  well  as  for  the 
cases  where  all  sensor  errors  and  no  sensor  errors  are 
present.  In  these  cases,  structural  modes  are  not  sensed. 

An  additional  simulation  is  used  to  illustrate  the  effects 
of  structural  modes.  The  flight  profile  and  sensor  para¬ 
meters  given  in  the  previous  chapter  are  used  in  the  gen¬ 
eration  of  these  error  time  histories.  The  analysis  and 
comparison  of  these  error  time  histories  are  presented 
in  three  major  sections «  effects  of  laser  gyro  errors, 
effects  of  accelerometer  errors,  and  effects  of  structural 
modes.  For  the  sections  on  sensor  error  effects,  a  general 
overview  of  the  error  time  history  obtained  with  all  errors 
present  will  be  presented,  followed  by  a  breakdown  of  the 
effects  of  the  individual  sensor  error  sources.  Following 
these  sections,  the  effects  of  the  structural  modes  will  be 
considered.  A  basis  for  the  evaluation  of  the  performance 


30 


of  the  INS,  when  structural  modes  are  not  sensed,  is  set 
at  -  3  ft/sec  velocity  error  and  -  1  nm/hr  position  error. 


Also,  in  analyzing  the  simualtion  results,  use  is 
made  of  the  error  state  equations  for  a  local  vertical 
platform  for  confirming  the  results. 

Effects  of  Laser  Gvro  Errors 

Figure  8  shows  the  time  history  obtained  for  the  laser 
gyro  with  all  errors  present.  The  roll  and  pitch  Euler 
angle  errors  and  the  east  and  north  velocity  errors  as 
well  as  the  longitude  and  latitude  errors  are  low  fre¬ 
quency  sinusoids  with  transients  due  to  the  maneuvers  super¬ 
imposed  upon  them.  The  yaw  angle  error  is  ramp-like  in 
nature  while  the  vertical  channel  errors,  altitude  and 
vetical  velocity,  are  relatively  unaffected  by  the  sensor 
errors.  The  reason  for  the  vertical  channel  being  indif¬ 
ferent  to  seuisor  errors  is  because  it  is  stabilized  with  a 
third  order  damping  system.  The  aircraft  maneuvers 
generally  result  in  Hie  Euler  angles  experiencing  step 
or  pulsed  transients.  The  local  level  velocity  errors 
are  effected  by  the  integration  of  the  effect  of  the  pitch 
and  roll  angle  errors  on  the  acceleration  vectors.  This 
effect  results  in  the  velocity  errors  experiencing  steps 
or  ramps.  These  same  effects  were  passed  on  into  the  posi¬ 
tion  errors,  but  since  the  position  errors  are  effected  by 
the  integration  of  the  velocity  errors,  the  effects  are 
more  benign.  The  position  errors  then  couple  back  to  effect 
the  roll  and  pitch  angle  errors.  Since  the  position  errors 


Figure  8c.  Simulation  Results  Obtained  with  no  Structural  Modes.  Ideal  Accelerometer 
and  all  Gyro  Errors  Present  (Attitude  Errors) 


have  little  effect  on  the  yaw  angle  error,  it  is  ramp-like 
in  nature.  One  of  the  most  prominent  effects  of  the  maneu¬ 
vers  is  that  of  the  loiter.  The  loiter  introduces  a  sinu¬ 
soidal  error  into  the  navigation  system  variables. 

The  time  histories  shown  in  figure  8  are  a  composite 
of  the  effects  of  each  individual  sensor  error.  Each  of 
these  error  sources  is  now  considered  separately  to  deter¬ 
mine  their  individual  effect  on  the  IKS  errors. 

Gvro  Bias  Error.  Figure  9  shows  the  plots  of  the  INS 
errors  for  simulations  with  only  bias  error  present,  all 
errors  present,  and  no  errors.  The  bias  error  plot  is  re¬ 
presented  by  the  square,  the  plot  of  all  errors  by  a  triangle* 
and  the  case  for  ideal  sensors  is  marked  with  a  square. 

The  presence  of  the  gyro  bias  error  results  in  the  Euler 
angles  initially  ramping  off  at  the  rate  of  the  bias  error 
(0.01°/hr),  The  yaw  angle  error  is  only  slightly  effected 
by  the  aircraft  maneuvers  and,  thus,  it  continues  to  ramp  off 
at  about  O.Oi°/hr.  Before  the  loiter  maneuver,  the  bias  error 
accounts  for  about  half  of  the  total  error,  but  during  the  loi¬ 
ter  the  bias  error  becomes  negligible  compared  to  the  other 
gyro  errors. 

The  snap  rolls  cause  a  step  in  the  bias  induced  pitch 
angle  error.  Also,  during  the  loiter  maneuver,  the  ampli¬ 
tude  of  the  oscillation  of  the  pitch  angle  decreases  more  ra¬ 
pidly  with  all  sensors  present  then  with  only  the  bias  error 
present.  This  indicates  that  the  bias  error  prolongs  the 
settling  of  errors  induced  by  the  maneuvers  of  the 
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aircraft.  The  bias  error  proved  to  be  the  primary  driver 
of  the  pitch  angle  error.  The  bias  error  induces  a  ramp 
of  0.01°/hr  in  the  roll  angle  error  that  is  only  slightly 
altered  by  any  aircraft  maneuvers  until  the  loiter  maneuver 
is  begun.  During  the  loiter,  as  in  the  case  of  the  pitch 
angle,  the  bias  induced  roll  angle  error  tends  to  slow  the 
settling  of  the  error  induced  in  the  roll  angle  by  the  man¬ 
euvers.  the  bias  error  accounts  for  a  large  part  of  the 
roll  angle  error  except  during  the  phases  of  the  flight 
where  the  snap  rolls  occur.  Here,  steps  are  introduced 
by  other  gyro  errors  and  this  causes  the  roll  angle  to  be¬ 
come  much  larger  when  all  gyro  errors  are  present  than  when 
the  bias  alone  is  present.  The  east  velocity  error  is  large¬ 
ly  influenced  by  the  integration  of  the  roll  angle's  effect 
on  the  vertical  acceleration.  This  gives  it  a  ramp-like 
nature  with  a  slope  of  approximately  (3.5  ft/sec )/hr.  The 
steps  are  introduced  into  the  roll  angle  when  all  sensor 
errors  are  present,  resulting  in  a  larger  slope  in  the 
east  velocity  error  with  all  gyro  errors  present  than  with 
only  the  bias  error  present.  This  indicates  that  the  bias 
error  is  not  the  main  contributer  to  the  east  velocity  error. 
The  effect  of  the  bias  error  on  the  north  velocity  error 
is  more  pronounced.  The  bias  error  contributes  a  major 
portion  of  the  north  velocity  error  due  to  gyro  errors. 

Also,  the  ramp  induced  in  the  north  velocity  error  by  the 
gyro  bias  is  only  slightly  effected  by  ’the  loiter  maneuver 
and  is  larger  than  the  error  due  to  all  gyro  errors. 
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The  position  errors,  latitude  and  longitude,  are 
mainly  driven  by  the  integral  of  the  velocity  errors. 
Therefore,  they  too,  exhibit  a  ramp-like  nature,  However, 
during  the  loiter  there  are  more  oscillations  in  the  position 
errors,  but  this  is  not  caused  by  the  gyro  errors  because 
these  oscillations  are  evident  even  when  ideal  sensors  are 
used.  As  is  expected,  since  position  errors  are  largely 
influenced  by  the  integration  of  the  velocity  errors,  the 
bias  error  has  a  similar  effect  on  the  position  errors  as 
it  does  on  the  velocity  errors. 

Gyro  Misalignment  Error.  The  second  gyro  source  ex¬ 
amined  is  the  misalignment  error.  Figure  10  shows  the 
.  results  of  the  simulations  made  with  misalignment  error 
only  plotted  against  simulation  results  for  the  ideal  gyro 
and  for  the  laser  gyro  with  all  errors  present.  The  most 
immediately  noticable  effect  of  the  misalignment  error 
occurs  in  the  yaw  angle  error  at  the  points  where  the  snap 
rolls  are  performed,  but  its  most  degrading  effect  occurs 
during  the  loiter.  The  misalignment  error  produces  a 
rather  steeply  sloping  yaw  angle  error  that  rises  until 
the  loiter  is  completed.  Appreciable  effects  from  the 
misalignment  error  are  also  found  in  the  pitch  and  roll 
angle  errors.  Steps  occur  in  these  errors  during  the  snap 
rolls  and  the  evasive  maneuvers. 

However,  the  step  in  the  pitch  angle  error  during  the 
second  snap  roll  is  in  the  opposite  direction  of  the  steps 
from  the  other  snap  rolls.  This  indicates  that  the  pitch 
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angle  is  sensitive  to  roll  direction,  while  the  roll  angle 
error  has  steps  at  each  snap  roll  that  are  in  the  same  di¬ 
rection.  The  misalignment  error  clearly  makes  up  for  a 
large  portion  of  all  of  the  Euler  angle  errors  and  during 
loiter  it  becomes  even  more  dominant.  The  step  errors  in 
the  roll  and  pitch  angle  errors  are  integrated  and  result 
in  ramp  errors  in  the  east  velocity  error.  This  ramp  error 
makes  up  for  most  of  the  difference  that  existed  between 
the  bias  error  induced  east  velocity  error  and  that  re¬ 
sulting  from  all  of  the  gyro  errors.  This  increased  ramp 
effect  is  enough  by  itself  to  push  the  east  velocity  error 
to  the  three  nautical  mile  error  value.  Therefore,  reduction 
of  the  misalignment  error  is  highly  desirable  for  improving 
the  INS  performance.  Because  of  the  large  yaw  error  and  the 
relatively  small  pitch  error,  the  north  velocity  error  is 
only  mildly  effected  by  the  misalignment  error  during  most 
of  the  flight.  (The  coupling  of  the  pitch  angle  to  vertical 
acceleration  tends  to  subtact  from  the  coupling  of  the  yaw 
error  with  the  east  acceleration).  However,  during  the 
loiter,  the  misalignment  error  causes  the  north  velocity 
error  to  oscillate  as  the  aircraft  turns  through  the  three 
turns.  Once  again,  since  the  position  errors  are  largely 
effected  by  the  integration  of  the  velocity  errors,  the 
effects  on  them  due  to  misalignment  error  are  similar  to 
those  of  the  velocity  errors. 

Gvro  Scale  Factor  Error.  A  third  laser  gyro  error 
source  considered  in  this  analysis  is  the  scale  factor 
error.  The  error  time  histories  obtained  for  it  are  shown 
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in  Figure  11.  The  scale  factor  has  less  overall  effect 
on  the  INS  errors  than  either  the  bias  or  the  misalign¬ 
ment  errors.  The  most  notable  effect  of  the  scale  factor 
error  is  the  steps  it  introduces  in  the  roll  angle  error 
when  the  aircraft  does  a  snap  roll.  Notice  that  the  se¬ 
cond  roll  creates  a  negative  step  error.  This  is  because 
the  second  roll  is  made  to  the  opposite  direction  of  the 
others.  The  value  of  these  steps  is  about  0.0019°.  Since 
the  second  roll  is  done  in  the  opposite  direction  of  the 
others,  the  effect  of  the  steps  on  the  east  velocity  error 
and  thus  the  longitude  error  caused  by  the  scale  factor 
is  small  compared  to  that  caused  by  the  misalignment  error. 

Gvro  Random  Walk  Error.  The  plots  obtained  to  illus¬ 
trate  the  effects  of  the  gyro  random  walk  error  on  the  INS 
performance  are  shown  in  Figure  12.  The  random  walk  error 
has  little  effect  on  the  INS  errors  when  compared  with  the 
effects  of  the  misalignment, bias,  and  scale  factor  errors. 
The  only  noticable  effect  the  random  walk  error  has  is  a 
step  error  in  pitch  angle  during  the  evasive  maneuver. 

This  effect,  however,  does  not  cause  any  appreciable  error 
in  the  velocity  and  position  parameters.  Therefore,  the 
random  walk  error  can  be  considered  negligible  when  using 
the  scale  factor,  bias,  and  misalignment  errors  specified. 

Summary .  Figure  13  shows  the  simualtion  results  for 
the  laser  gyro  with  all  errors  present  and  with  a  5^  min¬ 
ute  extension  of  level  flight  added  to  the  flight  profile. 
Adding  this  additional  flight  time  allowed  the  illustration 
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Case  with  ideal  sensors. 

Case  with  all  gyro  errors  present 
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with  all  gyro  errors  present 


of  the  errors  after  one  Schuler  period.  Examining  the  yaw 
angle  error,  it  is  seen  that  upon  completion  of  the  maneu¬ 
vers,  the  error  ramps  up  with  a  slope  equal  to  the  bias'  error 
just  as  it  did  at  the  outset.  The  pitch  and  roll  angle  ex¬ 
perience  oscillations  at  the  Schuler  frequency.  As  is  seen 
in  the  graphs,  the  amplitude  of  the  Schuler  frequency  os¬ 
cillation  of  the  roll  angle  is  smaller  than  the  largest 
value  of  step  errors  that  resulted  from  the  high  dynamic 
maneuvers.  Likewise,  the  pitch  angle  error  has  a  Schuler 
frequency  amplitude  that  is  smaller  than  the  largest  amp¬ 
litudes  of  the  oscillations  created  by  the  loiter.  The 
implications  of  this  are  that  the  strapped-down  system  de¬ 
finitely  experiences  degradation  in  performance  when  sub¬ 
jected  to  a  high  dynamic  environment. 

The  north  velocity  error  is  much  closer  to  meeting 
the  specifications  than  the  east  velocity  error.  Examin¬ 
ation  of  the  84  minute  simulation  shows  that  the  north 
velocity  error  reaches  a  peak  at  about  one  hour.  The  value 
at  this  point  is  -4  ft/sec.  This  results  in  a  latitude 
error  of  1  nra.  Thus  the  latitude  error  is  the  only  para¬ 
meter  that  reaches  the  specifications  and  it  is  quite  mar¬ 
ginal  . 

The  overall  effect  of  the  laser  gyro  errors  indicates 
that  improvements  must  be  made  before  a  performace  speci¬ 
fication  of  -  3  ft/sec  velocity  error  and^l  nm/hr  position 
error  can  be  achieved.  The  major  problems  are  in  east  velocity 
error,  and  thus  the  longitude  error.  These  error  result 
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from  the  large  roll  angle  errors  created  primarily  by  the 
misalignment  error,  while  the  bias  and  scale  factor  error 
would  have  had  about  twice  as  much  effect  because  the  sign 
of  the  step  errors  in  the  roll  angle  caused  by  the  scale 
factor  error  are  dependent  upon  the  direction  of  the  rolls. 
Conversely,  the  step  errors  induced  by  the  gyro  misalign¬ 
ment  error  added  regardless  of  the  roll  direction. 

Effects  of  Quartz-Flexure  Accelerometer  Errors 

A  similar  analysis  to  that  presented  on  the  laser 
gyro  errors  is  now  presented  for  the  accelerometer  errors. 
Figure  14  shows  the  error  time  histories  generated  for  the 
quarta-flexure  accelerometer  with  all  error  sources  present. 
Unlike  the  gyro  errors,  the  accelerometer  errors  produce  er¬ 
rors  in  the  velocity  and  thus  position  errors,  which  then 
couple  back  to  effect  the  velocity  errors.  However,  the 
loiter  maneuver  still  has  the  same  basic  effect  on  the 
accelerometer  error  induced  INS  errors.  That  is,  the  INS 
errors  oscillate  a6  the  aircraft  performs  the  loiter  maneu¬ 
ver. 

The  east  velocity  error  tends  to  be  relatively  stable 
except  for  step  errors  from  the  snap  rolls  and  evasive 
maneuver.  It  does  oscillate  during  the  loiter  maneuver, 
but  its  mean  remains  fairly  constant.  The  north  velocity 
error,  however,  tended  to  be  ramp-like  and  changed  directions 
at  the  point  of  the  start  of  the  loiter  maneuver.  The  posi¬ 
tion  errors  followed  the  general  nature  of  the  velocity 
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errors  except  that  the  latitude  error  did  not  change  direc¬ 
tions.  The  pitch  angle  error  resenbles  the  north  velocity 
while  the  roll  angle  is  of  the  same  form  as  the  east  velocity 
error.  The  yaw  angle  error  is  relatively  stable  except  for 
the  steps  induced  by  the  snap  rolls  and  evasive  maneuver. 

As  was  the  case  of  the  gyro  errors,  accelerometer  errors 
have  little  effect  on  the  vertical  channel.  The  effect 
of  the  individual  accelerometer  errors  on  INS  performance 
is  now  assessed. 

Accelerometer  Bias  Error.  The  bias  error  is  the  dom¬ 
inate  accelerometer  error  source.  Figure  15  shows  the  re¬ 
sults  of  the  simulation  obtained  with  the  bias  as  only 
accelerometer  error.  It  produces  a  ramp  in  the  east  vel¬ 
ocity  error  that  is  of  a  much  larger  magnitude  than  that 
for  the  case  where  all  accelerometer  errors  are  present. 

This  results  because  the  other  accelerometer  errors  create 
east  velocity  errors  that  are  of  the  opposite  sign.  Also, 
there  is  a  sign  change  in  the  east  velocity  error  at  the 
beginning  of  the  loiter.  The  longitude  error  follows  the 
same  general  characteristics  as  the  east  velocity  error 
except  that  it  does  not  change  signs.  Small  pulses  are 
visible  on  the  accelerometer  bias  error  induced  east  velo¬ 
city  error  at  the  points  of  the  evasive  maneuver  and  snap 
rolls.  The  pitch  and  roll  angles  both  ramp  off  before  os¬ 
cillating  as  the  loiter  maneuver  begins.  The  pitch  angle 
is  essentially  the  same  as  the  north  velocity  error.  The 
yaw  angle  error  with  bias  error  only,  basically  follows  the 


LRTITUOC  ERROR  LONGITUDE  ERROR  RLTITUOE  ERROR 


VELOCITY  ERROR  £flsr  VELOCITY  ERROR  ■  VERTICAL  VELOCITY 


same  wave  form  as  the  yaw  angle  error  with  no  accelero¬ 
meter  errors  present,  except  that  it  ramps  off  after  ex¬ 
periencing  a  step  error  at  the  outset  of  the  climb. 

Accelerometer  Misalignment  Error.  Figure  l6  shows 
the  effect  of  the  accelerometer  misalignment  error  on  the 
INS  errors.  Although  the  pulses  induced  at  the  maneuvers 
are  more  pronounced  and  the  overall  error  more  unsettled, 
the  east  velocity  error  that  results  from  accelerometer 
misalignment  has  the  same  general  form  as  that  for  the 
bias  error  only  in  the  opposite  direction.  The  magni¬ 
tude  is  about  six  times  smaller  than  that  due  to  bias,  but 
it  still  has  the  effect  of  offsetting  some  of  the  error 
caused  by  the  bias.  Also,  the  pulse  introduced  into  the 
east  velocity  error  by  the  evasive  maneuver  is  more  pro¬ 
nounced  with  the  misalignment  error.  The  longitude  error 
that  results  is  also  of  the  opposite  .sign  and  about  six 
times  smaller  than  that  achieved  with  the  accelerometer 
bias  error.  The  accelerometer  misalignment  had  little 
effect  on  the  north  velocity  error,  and  thus,  the  lati¬ 
tude  error  was  also  relatively  uneffected  by  misalign¬ 
ment.  The  misalignment  error  did  however,  have  an  appre¬ 
ciable  effect  on  the  roll  angle  error.  The  roll  angle 
experiences  steps  at  the  snap  rolls,  with  the  direction 
of  the  steps  being  dependent  on  the  direction  of  the  roll. 
This  indicates  that  if  the  second  roll  had  been  in  the  same 
direction  as  the  other  rolls,  the  error  in  roll  angle  pro¬ 
duced  by  the  misalignment  error  would  have  been  increased. 
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□  Case  with  accelerometer  misalignment  only. 
O  Case  with  ideal  sensors. 

A  Case  with  all  accelerometer  errors  present. 


Case  with  accelrometer  misalignment  error  only 
Case  with  ideal  sensors. 

Case  with  all  accelerometer  errors  present. 


The  evasive  maneuver  introduces  a  pulse  that  is  about  half 
the  magnitude  of  the  pulse  produced  when  all  accelerometer 
errors  are  present.  Pitch  and  yaw  angles  show  negligible 
effects  from  the  misalignment  error. 

Accelerometer  Scale  Factor  £rror.  The  effects  of  the 
accelerometer  scale  factor  error  on  the  INS  parameter  is 
illustrated  in  Figure  17.  The  north  velocity  error  is  only 
slightly  effected  by  the  accelerometer  by  the  accelerometer 
scale  facotr  error.  Likewise  the  scale  factor  error  produces 
negligible  error  in  latitude.  However,  in  regards  to  east 
velocity  error,  the  scale  factor  error  accounts  for  a 
large  portion  of  this  error  during  the  period  before  the 
second  snap  roll  and  during  the  loiter  maneuver.  The 
reason  for  the  decreased  effect  of  the  scale  factor  error 
between  the  second  snap  roll  and  the  loiter  results  be¬ 
cause  the  error  in  the  east  velocity  due  to  the  scale 
factor  error  is  in  the  opposite  direction  from  a  similar 
step  error  caused  by  some  other  source.  The  longitude 
error  exhibited  similar  characteristics  as  the  east  velo¬ 
city  error.  The  large  spike  that  occurs  at  the  evasive 
maneuver  which  appeared  as  if  it  may  be  caused  by  the  scale 
factor  error,  is  much  clearer  in  the  longitude  error.  Since 
it  is  present  for  the  case  of  an  ideal  accelerometer,  it  must 
be  caused  by  some  source  other  than  accelerometer  errors. 

The  Euler  angles  are  all  relatively  uneffected  by  the  scale 
factor  error. 
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Case  with  all  accelerometer  errors  present. 


Other  Accelerometer  Errors.  The  other  two  error  sources 
modeled  for  the  quartz-flexure  accelerometer  are  the  non¬ 
linearity  and  cross-coupling  error.  Figures  18  and  19  show 
the  plots  obtained  for  each  of  these  cases.  The  only  no- 
ticable  effect  either  error  source  has  on  the  INS  errors  is 
the  effect  of  the  nonlinearity  error  on  the  east  velocity 
error  during  the  period  from  the  last  snap  roll  until  the 
loiter.  During  this  period,  the  east  velocity  error  due 
to  accelerometer  non-linearity  subtacts  from  the  velocity 
error  that  results  when  no  sensor  errors  are  present.  This 
creates  a  small  reduction  in  the  overall  error  in  east 
velocity  and  thus  the  longitude  error  during  this  phase 
of  the  flight. 

Summary .  The  overall  errors  in  the  INS  parameters 
caused  by  the  accelerometer  errors  are  somewhat  smaller 
than  those  resulting  from  laser  gyro  errors.  The  bias  and 
misalignment  error  dominate  the  effect  of  the  accelerometer 
errors.  A  reduction  in  the  accelerometer  errors  should  re¬ 
sult  in  a  direct  reduction  of  the  velocity  errors  since  the 
accelerometer  errors  are  coupled  into  the  velocities  in  an 
additive  fashion.  Since  the  position  and  attitude  errors 
are  driven  by  the  velocity  errors,  reduction  in  the  velocity 
errors  should  result  in  reductions  in  both  position  and 
attitude  errors. 

As  is  done  for  the  case  of  laser  gyro  errors,  an  84 
minute  simulation  is  made  to  illustrate  the  effects  of  the 
accelerometer  errors  over  one  Schuler  period.  Figure  20 
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These  are  so  close  in  magnitude  as  to  be  indistinguishable. 


Simulation  Results  Obtained  Over  one  Schuler  Period  for  Case 
Accelerometer  Errors 


1 


shows  the  error  time  histories  that  are  obtained  from 
this  simulation.  Once  again,  as  is  true  for  the  case  of 
laser  gyro  errors,  the  major  portion  of  the  errors  is  a  re¬ 
sult  of  the  effects  of  the  high  dynamic  maneuvers  performed 
during  the  first  30  minutes  of  flight.  However,  The  values 
of  the  INS  errors  due  to  accelerometer  errors. are  much 
smaller  than  those  caused  by  the  laser  gyro  errors.  Also, 
both  east  and  north  velocity  errors  are  within  the  specifi¬ 
cation  as  is  the  position  errors.  This  indicates  that  re¬ 
ducing  the  laser  gyro  errors  should  be  more  beneficial 
than  reducing  accelerometer  errors. 

Effect  of  Structural  Modes 

The  simulations  discussed  in  the  previous  sections 
do  not  account  for  structural  mode  effects.  This  section 
presents  simulations  obtained  with  structural  modes  sensed 
and  provides  an  analysis  of  their  effects  on  the  INS  errors. 

Figure  21  shows  the  time  histories  of  the  INS  errors 
with  structural  modes  sensed  plotted  against  those  with 
no  structural  modes  sensed.  In  both  cases  the  sensors 
are  assumed  ideal. 

The  structural  modes  produce  a  definite  degradation 
in  the  overall  INS  performance.  The  roll  angle  is  effected 
less  than  the  other  attitude  errors.  With  structural  modes 
present,  the  roll  and  evasive  maneuvers  produce  some  large 
pulse  transients  into  the  roll  angle.  The  one  from  the  eva¬ 
sive  maneuver  is  nearly  twice  as  large  as  those  from  the 
rolls.  Also,  a  sinusoid  results  during  the  loiter. 
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Part  of  the  cause  of  this  sinusoid  is  the  coupling  of 
the  latitude  error  into  the  roll  angle  error.  The 
coupling  of  the  other  attitude  errors  also  adds  to  the 
sinusoid. 

The  presence  of  the  structural  inodes  causes  the  pitch 
angle  to  be  quite  sensitive  to  the  aircraft  maneuvers. 

The  pitch  angle  exhibits  a  ramp  during  the  initial  des¬ 
cent.  At  each  of  the  snap  rolls,  the  pitch  angle  error 
experiences  steps  on  the  order  of  0.03°.  Even  though  one 
of  the  rolls  is  in  the  opposite  direction  of  the  other  two, 
the  steps  caused  by  the  snap  rolls  have  an  additive  effect. 
The  evasive  maneuver  produces  a  negative  step  error  that  is 
about  twice  as  large  as  that  of  the  snap  roll.  The  pitch 
angle  error  is  also  sensitive  to  the  weapon  delivery.  A 
pulse  transient  results  from  this  maneuver.  Although  much 
smaller  than  the  one  at  the  weapon  delivery,  the  ascent 
produces  a  transient  in  the  pitch  angle  error,  also.  These 
transients  reflect  the  shape  of  the  longitudinal  structural 
modes.  During  the  periods  of  flight  where  the  pitch  angle 
error  is  not  experiencing  steps  or  transients,  it  exhibited 
a  ramp-like  nature.  Also,  as  in  the  roll  angle  error,  the 
latitude  error  couples  back  into  the  pitch  angle  error  during 
the  loiter. 

The  yaw  angle  also  has  significant  errors  induced  into 
it  by  the  aircraft  maneuvers  when  the  structural  modes  are 
sensed.  However,  the  yaw  angle  error  that  results  is  nega¬ 
tive.  This  should  have  the  effect  of  reducing  the  yaw  angle 
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error  when  the  laser  gyro  misalignment  error  is  present. 

The  roll  angle  error  couples  into  the  yaw  angle  error  re¬ 
sulting  in  large  transients  at  each  snap  roll.  The  evasive 
maneuver  results  in  a  large  negative  step  error.  After  the 
final  snap  roll  has  taken  place,  the  yaw  angle  becomes  quite 
benign.  There  is  only  a  small  sinusoid  during  the  loiter, 
therefore,  there  must  be  a  small  amount  of  coupling  from 
the  latitude  errors. 

One  of  the  most  significant  effects  of  the  structural 
modes  is  in  the  north  velocity  error.  The  north  velocity 
reaches  a  value  of  -18  ft/sec.  The  steps  that  result  in  the 
attitude  angle  errors  during  the  aircraft  maneuvers  inte¬ 
grate  into  the  north  velocity  error  causing  the  ramp-like 
behavior.  As  a  result  of  this  velocity  error,  the  lati¬ 
tude  error  is  also  driven  to  exhibit  a  ramp-like  nature  and 
is  quite  large  (3.6  nm)  by  the  end  of  the  30  minute  flight. 

The  east  velocity  error  is  less  effected  by  the  struc¬ 
tural  modes.  Until  the  loiter,  the  only  significant  effect 
of  the  structural  modes  occurs  when  the  aircraft  performs 
the  evasive  maneuver.  Here  a  transient  occurs  that  has  a 
peak  value  of  about  0.5  ft/sec,  However,  during  the  loiter, 
the  east  velocity  error  begins  to  ramp  off,  settling  again 
after  completion  of  the  loiter.  The  peak  error  during  this 
period  is  less  than  2  ft/sec.  The  longitude  error  then 
follows  the  same  characteristics  as  the  east  velocity 
error,  but  the  transients  are  more  benign  due  to  the  inte¬ 
gration. 
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The  structural  inodes  also  have  some  interesting  effects 
on  the  vertical  channel  of  the  INS.  The  vertical  velocity 
error  experiences  much  larger  transients  from  the  maneuvers 
when  structural  modes  are  present.  The  largest  ones  are  at 
the  snap  rolls  and  evasive  maneuver.  Also,  during  the  loiter 
the  amplitude  of  the  transients  are  quite  large  when  struc¬ 
tural  modes  are  present  as  opposed  to  when  they  are  not.  The 
vetical  velocity  error  does  however  remain  fairly  small, 
with  the  largest  transient  being  smaller  than  i  ft/sec. 

The  transients  of  the  attitude  error  never  exceeded  20 
feet  while  the  mean  remains  nearly  zero. 

Summary .  The  overall  effects  of  the  structural  modes 
on  the  INS  performance  are  not  good.  The  most  significant 
effect  of  the  structural  modes  occurs  in  the  north  velocity 
error.  Here  the  structural  modes  cause  an  error  of  nearly 
20  ft/sec.  This  then  results  in  a  position  error  of  close 
to  4nm.  These  errors  are  much  larger  than  the  position 
and  velocity  errors  caused  by  sensor  errors.  Thus,  no 
matter  how  much  improvement  is  made  on  the  effects  of  the 
sensor  errors,  the  sensors  must  be  located  such  that  the 
structural  mode  effects  are  minimized  before  a  desirable 
level  of  performance  can  be  achieved  from  the  strapped- 
down  INS  in  a  high  dynamic  environment. 
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Conclusions 

In  this  chapter,  the  effects  of  the  sensor  errors  and 
the  structural  modes  on  INS  performance  have  been  assessed. 
The  analysis  reveals  that  the  high  dynamic  environment 
results  in  significant  degradation  in  the  INS  performance. 
The  INS  errors  in  local  level  position  and  velocity  caused 
by  laser  gyroerrors  are  found  not  to  meet  the  specifications 
of  *  3  ft/sec  velocity  error,  while  those  caused  by  the  accel 
erometer  errors  are  within  the  specifications.  This  indi¬ 
cates  that  improvements  must  be  made  in  the  laser  gyros 
and/or  the  methods  used  to  sense  the  rotation  rates.  A 
change  of  gyro  configuration,  for  example,  could  result  in 
less  effect  from  the  coupling  of  the  roll  angle  error  into 
the  east  velocity  error. 

From  examining  the  effect  that  structural  modes  have 
on  the  INS  errors  it  is  found  that  they  produce  additional 
degradation  to  that  of  the  sensors.  The  most  significant 
effect  of  the  structural  modes  is  the  large  pitch  angle 
errors  that  coupled  into  the  north  velocity  causing  a  se¬ 
vere  error  there.  Also,  since  the  errors  caused  by  the 
structural  modes  basically  added  to  those  caused  by  the 
sensor  errors,  and  since  the  gyro  induced  errors  were  lar¬ 
ger,  the  location  of  the  gyro  is  more  critical  than  that  of 
the  accelerometers. 
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The  sensitivity  of  the  INS  errors  to  changes  in  the 


values  of  the  sensor  errors  are  determin  ed  in  this  chap¬ 
ter.  A  simulation  is  made  with  each  sensor  error  reduced 
to  50%  of  the  value  used  for  that  parameter  in  the  error 
analysis  presented  in  Chapter  3*  The  results  are  then 
plotted  against  the  siraualtion  made  with  all  sensor  errors 
present  with  the  parameter  values  as  defined  in  Chapter  two. 
These  plots  are  then  compared  to  determine  the  sensitivity 
of  the  INS  errors  to  determine  the  sensitivity  of  the  INS 
errors  to  each  sensor  error. 

The  results  of  this  analysis  pertaining  to  the  laser 
gyro  errors  are  presented  first.  Following  this,  the  same 
type  of  analysis  is  presented  for  the  quartz-flexure  accel¬ 
erometer. 

Sensitivity  to  Laser  Gyro  Errors 

The  plots  in  Figure  22  show  the  effects  of  reducing  the 
laser  gyro  bias  error  by  $0 %.  The  curve  labeled  with  a  cir¬ 
cle  represents  simualtion  results  with  all  laser  gyro  errors 
present,  while  the  curve  labeled  with  a  triangle  represents 
the  results  with  no  laser  gyro  errors.  The  third  curve, 
labeled  with  a  square,  is  the  results  obtained  with  all 
laser  gyro  errors  present,  but  the  bias  error  at  a  value  of 
50#  less  than  the  one  used  in  previous  simulations  (0.005°/hr). 


Figure  22a.  Comparison  of  Position  Errors  for  Cases  of  All  Laser  Gyro  Errors 
Full  Value,  Gyro  Bias  Error  at  5°^  Value,  and  all  Ideal  Sensors 


£2 

VT\ 

H 

+*1—4  • 

rf  3  ro 
o 

o  +>  w 
U  ri  C 
U  o> 
0)  n  to 

M  OH 
ri  ri 

•H  t-i  <1) 

jo  a> 

•H 

>»o 

i— i  i— i 

C  SH 

o  fcD  ri 

jo  jo  x: 
+*+>+> 
•rl  *r|  *r4 

•s  5  ? 

a,  o>  a> 
rot o  w 
ri  ri  ri 
ouo 


□O0 


Recalling  the  analysis  of  chapter  three  in  regards  to 
the  bias  error,  it  should  be  expected  that  a  change  in  the 
value  of  bias  error  should  have  little  effect  on  the  yaw  an¬ 
gle  error.  This,  indeed,  is  the  case.  The  50%  reduction 
in  bias  error  resulted  in  a  nearly  negligible  reduction  in 
yaw  angle  error.  However,  the  step  errors  in  the  roll  an¬ 
gle  are  slightly  reduced  and  since  the  steps  are  additive, 
additional  improvement  over  the  results  obtained  with  the 
0.01°/hr  bias  error  are  obtained  at  each  snap  roll  and  the 
evasive  maneuver.  This  results  in  about  a  10%  improvement 
in  the  roll  angle  error  at  the  point  where  the  loiter  begins. 
During  the  loiter,  the  effects  of  misalignment  became  domi¬ 
nant  and  negligible  improvement  is  gained  by  a  5°%  reduction 
in  bias  error.  The  pitch  angle  error  shows  the  most  improve¬ 
ment  of  all  the  Euler  angles  with  the  reduction  in  bias  error. 
It  too  however,  is  not  greatly  improved,  except  during  climb. 
During  this  phase  of  the  flight  the  5 0%  reduction  in  bias 
error  results  in  a  50%  reduction  in  pitch  angle  error. 

This  is  a  good  illustration  of  the  fact  that  the  effects 
of  the  bias  error  are  less  sensitive  to  high  dynamic  maneuvers 
than  those  of  the  other  gyro  errors. 

With  the  yaw  and  roll  angle  errors  being  only  slightly 
reduced,  the  east  velocity  error  has  a  small  reduction  in 
the  value  of  its  slope.  By  the  end  of  the  30  minute  flight, 
this  has  resulted  in  a  reduction  of  the  error  by  0.6  ft/sec. 
This  indicates  that  the  specification  ofi3  ft/sec  position 
error  could  not  be  reached  by  simply  reducing  the  bias  error. 
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This  small  improvement  in  east  velocity  error  results  in  an 
improvement  in  the  longitude  error  of  0,0018°,  which  is  not 
enough  to  bring  it  to  within  the  specification  *  lnm/hr. 

As  is  found;  in  the  analysis  of  chapter  three,  the 
bias  error  results  in  most  of  the  north  velocity  error  when 
all  errors  of  the  lase  gyro  are  present.  This  indicates  that 
the  north  velocity  error  should  he  quite  sensitive  to  bias 
error  changes.  As  seen  in  Figure  22,  this  is  indeed  the 
case.  The  slope  of  the  north  velocity  error  is  reduced  such 
that,  by  the  end  of  the  30  minute  flight,  it  is  almost  half 
of  the  value  that  is  achieved  with  a  0.0l°/hr  bias  error. 

This  brought  the  value  of  the  errordown  to  abcrut  -0.6  ft/sec 
after  30  mintes  of  flight,  therefore  the  3  ft/sec  velocity 
error  specification  should  easily  be  reached. 

This  improvement  in  the  north  velocity  error  is  passed 
on  to  the  latitude  error.  After  30  minutes  of  flight,  the 
latitude  error  is  -0.0012  degrees  of  0.072  nm.  This.too, 
should  easily  remain  within  the  specification  of  1  nm/hr 
position  error. 

Since  the  vertical  channel  proved  to  be  relatively  un¬ 
effected  by  sensor  errors,  improvements  in  the  sensor  errors 
have  little  or  no  effect  on  the  altitude  and  altitude  rate 
errors . 


From  the  analysis  of  chapter  three,  the  misalignment 


error  is  found  to  effect  the  yaw  and  roll  angle  errors  and 


the  east  velocity  error  more  than  the  other  sensor  errors. 
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Figure  23  compares  the  INS  time  histories  achieved  with  the 
misalignment  error  at  5x10“-’  radians.  The  yaw  angle  error  is 
reduced  by  nearly  50, i,  while  the  pitch  angle  is  hardly  effect¬ 
ed.  The  roll  angle  experiences  improvements  at  the  points 
where  the  misalignment  error  intorduces  steps  into  it. 

This  occurs  at  the  snap  rolls  and  the  evasive  maneuver.  The 
improvement  that  results  is  about  n50;'-  reduction  in  the 
amount  of  step  error  caused  by  tho  misalignment  error. 

This  results  in  the  roll  angle  error  being  reduced  by  this 
amount  at  each  of  the  snap  rolls,  and  while,  at  the  evasive 
maneuver,  the  roll  angle  error  is  cut  in  half.  However,  dur¬ 
ing  tho  loiter,  this  improvement  is  reduced  such  that  by 
the  end  of  tho  flight  there  is  no  more  improvement.  Once 
again,  this  illustrates  the  fact  that  the  misalignment 
error  is  much  more  sensitive  to  aircraft  maneuvers  than  the 
bias  error. 

Since  the  yaw  and  roll  angle  errors  couple  into  the 
east  volocity,  by  tho  end  of  the  30  minute  flight,  the 
slope  change  has  resulted  in  tho  velocity  orror  being 
reduced  by  about  2  ft/soc.  This  puts  it  just  within  the 
specification  that  is  used  for  this  study.  The  east  velo¬ 
city  improvement  integrates  into  an  improvement  in  the  longi¬ 
tude  error.  This  improvement  is  ubout  0.00,  degrees  or 
0.3  nautical  miles.  This  should  bo  enough  of  an  improve¬ 
ment  ot  keep  tho  latitudo  error  loss  than  1  nm/hr. 

Since  the  misalignment  error  is  found  to  have  littlo 
offoct  on  north  volocity  and  latitudo  orrorn,  littlo  sensi¬ 
tivity  to  gyro  misalignment  error  is  expected.  The  simulation 
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results  obtained  with  the  misalignment  error  reduced  50$ 
verify  this  idea.  The  improvement  in  north  velocity  error 
is  less  than  0.25  ft/sec,  while  the  latitude  error  is  only 
reduced  by  10“  degrees  (6x10"-'  nautical  miles). 

Laser  Gyro  Scale  Factor  Error 

Unlike  the  bias  and  misalignment  errors,  the  scale 

factor  error  is  not  a  dominant  force  in  any  of  the  INS  errors. 
However,  it  does  account  for  about  V3  of  the  step  error 
size  in  the  roll  angle  error  and  has  about  the  same  effect 
as  the  bias  error  or  yaw  angle  error.  These  effects  couple 
into  local  level  velocity  and  position  errors,  with  more 
effect  being  on  the  east  than  the  north  velocity.  As  can 
be  seen  in  Figure  24,  when  the  scale  factor  error  is  re¬ 
duced  by  50$,  improvements  are  no ti cable  in  both  the  yaw 
and  roll  angles  and  these  improvements  help  to  reduce  the 
local  level  velocity  and  position  errors. 

The  improvements  in  the  yaw  angle  error  are  small  when 
compared  with  the  overall  yaw  angle  error,  but  the  improve¬ 
ments  in  roll  angle  error  are  substantial  during  the  snap 
rolls  and  evasive  maneuvers.  Since  both  the  roll  and  yaw 
angle  errors  couple  into  the  east  vleocity  error,  the  re¬ 
duction  of  scale  factor  error  results  in  almost  2  ft/sec 
improvement  in  east  velocity  error  by  the  end  of  the 
30  minute  flight.  This  translates  into  a  longitude  error 
of  0.004  degrees  (0.24  nautical  miles).. 

The  improvements  in  the  north  velocity  and  latitude 
errors  are  not  as  substantial  since  the  roll  angle  does 
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NORTH  VELOCITY  ERR® 


not  couple  into  these  errors.  The  slight  decrease  in  yaw 
angle  error,  however,  does  reduce  the  slope  of  the  north 
velocity  error  resulting  in  a  0.25  ft/sec  reduction  of 
its  value.  This  improvement  results  in  a  2.5x10  de¬ 
gree  reduction  in  the  latitude  error. 

Laser  Gvro  Random  Walk  Error 

As  should  be  expected  after  the  analysis  of  Chapter 
three,  improvement  in  the  random  walk  error  does  little  to 
improve  the  INS  errors.  Looking  at  the  graphs  in  Figures 
25,  the  curves  plotted  with  all  sensor  errors  at  the  values 
specified  in  Table  2  are  not  distinguishable  from  those  ob¬ 
tained  when  the  random  walk  was  at  50#  of  this  value  (  0.0025 
degrees/hr) . 

Summary  of  Sensitivity  to  Laser  Gyro  Errors.  Having 
determined  the  sensitivity  of  each  individual  laser  gyro 
error,  an  84  minute  simualtion  is  made  with  all  gyro  errors 
reduced  50#.  Figure  26  shows  the  results  that  are  obtained. 
As  should  be  expected,  since  it  has  third  order  damping, 
the  vertical  channel  is  relatively  insensitive  to  changes 
in  the  gyro  error  paramenters.  However,  all  other  INS  errors 
show  about  a  reduction  as  a  result  of  reducing  the  sen¬ 
sor  errors  by  50#.  Table  4  has  a  listing  of  the  amount  by 
which  each  INS  error  is  reduced  when  a  given  gyro  error  is 
reduced  by  50#.  As  can  be  seen,  the  results  for  the  case 
of  reducing  all  gyro  errors  is  not  the  summation  of  the 
case  of  reducing  all  gyro  errors  is  not  the  summation  of 
the  cases  of  reducing  the  gyro  errors  individually.  This 
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Case  with  Qyro  errors  at  nu 
Case  with  all  sensors  ideal. 


reduced  50%* 
at  full  value 


all  reduced  50% 
all  at  full  val 


TABLE  IV 


to  be  expected  since  the  Euler  angle  errors,  which  are 
directly  related  to  gyro  errors,  create  errors  in  the  ac*- 
celerometer  measurements  which  do  not  necessarily  have  the 
same  sign.  This  indicates  that  the  effects  of  the  gyro 
errors  can  be  overlapping  or  offsetting  to  some  extent 
depending  upon  the  direction  of  the  acceleration  along 
a  given  axis. 

The  simulations  presented  in  this  chapter  clearly  sjjow 
that  the  INS  errors  are  quite  sensitive  to  the  laser  gyro 
errors  when  structural  modes  are  not  sensed.  Table  4  shows 
the  results  of  the  sensitivity  tests.  Examining  the  figures 
presented,  it  is  seen  that  the  east  velocity  and  thus  the 
longitude  would  be  greatly  reduced  by  a  decrease  in  either 
the  scale  factor  error  or  the  misalignment  error,  while  a 
reduction  in  north  velocity  and  latitude  requires  an  im¬ 
provement  in  the  bias  error.  Since  the  east  velocity  and 
longitude  errors  are  found  to  be  the  worst,  the  misalign¬ 
ment  and  scale  factor  errors  should  be  given  more  attention 
than  the  bias  error. 

Sensitivity  to  Accelerometer  Errors 

A  similar  analysis  to  that  presented  on  the  sensiti¬ 
vity  of  INS  errors  to  laser  gyro  errors  is  now  presented, 
for  the  accelerometer. 

Accelerometer  Bias  Error.  The  INS  error  time  histories 
obtained  with  the  accelerometer  bias  error  reduced  to  50# 
are  shown  in  Figure  27*  As  is  shown  in  Chapter  three,  the 
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accelerometer  errors  effect  the  velocity  errors  directly. 

When  the  bias  error  is  reduced  50ft,  the  north  velocity  error, 
after  the  30  minute  flight,  is  reduced  by  only  a  small  amount. 
However,  just  before  the  loiter  maneuver  is  begun  the  smaller 
bias  error  provides  a  noticable  decrease  in  the  north  velo¬ 
city  error,  (about  27ft),  which  is  passed  on  to  the  latitude 
error  in  the  form  of  a  ramp  of  reduced  slope.  This  results 
in  a  reduction  in  the  latitude  error  of  about  the  same  amount. 
About  the  same  amount  of  improvement  in  the  pitch  angle  error 
should  be  expected  and  indeed  this  is  the  case. 

Since  the  bias  error  causes  a  large  east  velocity  and 
longitude  error  that  is  offset  by  the  effects  of  the  mis¬ 
alignment  and  scale  factor  error,  reducing  the  bias  error 
by  50ft  results  in  a  large  error  (10  times  that  with  all  errors 
at  full  value  for  the  east  velocity  and  8  times  that  with  all 
errors  for  longitude)  in  these  two  INS  paramenters.  However, 
this  error  is  the  negative  of  that  produced  when  only  the 
bias  error  is  present  since  it  is  due  to  the  other  acceler¬ 
ometer  gyro  errors.  This  indicates  that  if  the  accelerometer 
errors  are  to  be  reduced,  the  reduction  in  the  bias  error  must 
be  a  comparable  reduction  in  the  other  acceleometer  errors. 

Accelerometer  Misalignment  Error.  The  misalignment 
error  has  the  same  type  of  effect  on  the  INS  errors  as  the 
bias  error  except  that  the  effect  on  the  east  velocity  and 
longitude  are  opposite  in  sign  and  about  l/6  of  the  magni¬ 
tude  as  those  of  the  bias  error.  As  seen  in  Figure  28, 
this  results  in  the  north  velocity  and  latitude  errors  having 
about  the  same  sensitivity  to  misalignment  error  as  they  do 
to  bias  error,  2?ft.  The  east  velicity  and  longitude  errors  are 
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Comparison  of  Velocity  Errors  for  Cases  of  all  Accelerometer  Errors 
at  Full  Value,  Accelerometer  Bias  Reduced  and  all  Ideal  Sensors 


Case  with  ideal 
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RN6LE  ERROR 


effected  about  the  same  as  when  the  bias  is  reduced,  except 
that  the  errors  are  of  opposite  sign. 


Other  Accelerometer  Errors.  The  other  accelerometer 
errors,  scale  factor,  cross-coupling,  and  nonlinearity,  have 
little  effect  on  INS  errors  and  therefore  the  sensitivity  of 
the  INS  errors  to  these  errors  is  practically  zero.  The 
results  obtained  for  these  are  shown  in  Figures  29,  30,  and  31* 

Summary  of  Sensitivity  to  Accelerometer  Error.  The  most 
notable  effect  of  the  sensitivity  of  the  INS  errors  to  the 
accelerometer  errors  is  in  the  way  the  bias  and  misalignment 
errors  offset  each  others  effects  on  the  east  velocity  and 
longitude  errors  increasing  10  times  and  8  times  respectively 
when  either  the  bias  or  misalignment  error  is  reduced.  This 
indicates  that  it  is  desirable  to  make  comparable  changes  in 
both  accelerometer  errors  at  the  same  time  rather  than  just  one 
or  the  other. 

The  north  velocity  and  latitude  errors  are  equally  sen¬ 
sitive  to  changes  in  bias  and  misalignment  error.  However, 
the  effective  change  is  only  2?;i  for  each.  (Figure  32  shows 
a  comparison  of  the  results  obtained  when  all  accelerometer 
errors  were  reduced  and  those  obtained  with  the  original 
values  for  an  84  minute  simualtion).  Meanwhile,  the  east 
velocity  and  longitude  errors  are  increased  when  the  accel¬ 
erometer  bias  and  misalignment  errors  are  reduced  separately, 
because  these  INS  parameters  are  dependent  upon  the  difference 
in  these  two  errors  rather  than  the  sum. 
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NORTH  VELOCITY  ERROR 
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NORTH  VELOCITY  ERROR  EflST  v&-OCITY  ERROR  VERTICAL  VELOCITY  ERROR 
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Case  with  accelerometer  errors  at  full  value. 

Case  with  ideal  sensors. 

*These  are  so  close  in  magnitude  as  to  be  essentially  the  same  curves. 


Case  with  accelerometer  cross-coupling  error  reduced  50£. 

Case  with  accelerometer  errors  at  full  value. 

Case  with  ideal  sensors. 

*  These  are  so  close  in  magnitude  as  to  he  essentially  the  same  curves 


Case  with  accelerometer  nonlinearity  error  reduced  50#. 

Case  with  accelerometer  errors  at  full  value. 

Case  with  ideal  sensors. 

♦These  are  so  close  in  magnitude  as  to  be  essentially  the  same  curve 
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Figure  31c.  Comparison  of  Attitude  Errors  for  Cases  of  all  Accelerometer 
Errors  at  Full  Value,  Accelerometer  Nonlinearity  Reduced  50# 
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Table  5  shows  a  comparison  of  the  results  obtained 
in  the  sensitivity  analysis  on  the  accelerometer.  Note 
that  when  the  bias  and  misalignment  error  are  reduced  sep¬ 
arately!  the  east  velocity  and  longitude  errors  are  in-  .  ^  ■> 
creased.  This  increase  disappears  when  both  error  sources 
are  reduced  since  the  east  velocity  error  and  thus,  longi¬ 
tude  are  dependent  upon  the  difference  in  the  two  gyro  errors. 

Conclusions  on  Sensitivity  Analysis 

After  examining  the  effect  on  INS  performance  of  both 
l&sbr  gyro  and  accelerometer  errors,  it  is  found  that  the 
laser  gyro  errors  create  more  degradation.  The  findings 
show  that  improvements  in  gyro  errors  are  necessary  to  a- 
chieve  a  performace  specification  of  ^  lnn\/hr  position  error 
and  ^  3  ft/sec  velocity  error  when  no  structural  modes  are 
sensed . 

The  sensitivity  analysis  just  discussed  indicates 
that  this  performance  specification  is  achievable  by  reduc¬ 
ing  the  sensor  errors  by  50^.  Figure  33  shows  the  INS  error 
time  histories  for  the  case  using  the  sensor  error  values 
stated  in  Chapter  two,  Tables  1  and  2,  and  that  using 
50#  of  these  values.  This  shows  that  with  the  50#  values 
for  the  sensor  errors,  the  performance  specification  is  met. 
This  specification  could  also  be  met  at  some  lower  value, 
or  a  differnt  combination  of  lower  values,  of  the  sensor 
errors . 

The  implications  of  this  analysis  have  been  that  to 
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TABLE  V 


LATITUDE  ERROR  LONGITUDE  ERROR  ALTITUDE  ERROR 


bring  strapped-down  INS  performance  in  a  high  dynamic 
environment  down  to  a  competitive  level  with  gimbaled 
systems,  improvements  in  the  quality  of  gyros  used  on  the 
platforms  is  a  necessary  step. 
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Conclusion;;  nnd  Recommendations 

This  section  consists  of  two  parts.  The  first  part 
is  the  conclusions  that  have  been  drawn  from  this  research 
effort  and  the  sencond  is  recommendations  for  further  study 
in  this  research  area. 

Conclusions 

I 

Several  conclusions  can  be  drawn  from  this  research 
effort.  First  of  all*  it  has  been  shown  that  the  laser  gyro 
offers  better  performance  than  the  mechanical  two-dcgree-of- 
freedom  gyro.  This  is  determined  by  comparing  results  ob¬ 
tained  using  presently  achievable  performance  specifications 
for  each  of  these  classes  of  instruments.  The  INS  perfor¬ 
mance  obtained  for  each  is  then  compared  against  a  perfor- 
amnce  specification  of  -  1  nm/hr  position  error  and  -3  ft/sec 
velocity  error  when  no  structural  modes  are  sensed.  Due  to 
its  inherently  larger  misalignment  error,  the  two-degree- 
of-freedom  gyro  was  found  to  result  in  a  much  larger  east  ve¬ 
locity  and,  thus,  longitude  error,  than  the  laser  gyro. 

Another  conclusion  that  is  drawn  from  this  work  is  in 
regards  to  the  study  on  the  individual  sensor  error  and  struc¬ 
tural  mode  effects  on  INS  performance.  In  this  portion  of 
this  work,  it  is  found  that  the  INd  errors  are  degraded  more 
so  by  the  laser  gyro  errors  than  the  accelerometer  errors. 


Also,  it  is  found  that  the  structural  modes  create  substan¬ 
tial  degradation,  especially  to  the  north  velocity  and  lati¬ 
tude  errors,  of  the  INS  performance.  This  analysis  is  made 
based  on  the  assumption  that  the  sensor  are  in  an  orthogonal 
triad  at  the  center-of-gravity  of  the  aircraft.  This  result, 
coupled  with  the  laser  gyro  errors  being  more  of  a  factor  in 
INS  performance  than  accelerometer  errors,  indicates  that 
strong  consideration  of  the  location  of  the  gyro  should  be 
made. 

Still  another  important  point  brought  out  in  this  re¬ 
search  is  the  sensitivity  of  the  INS  errors  to  sensor  errors. 
Since  the  effects  of  the  accelerometer  errors  on  INS  per¬ 
formance  are  practically  negligible  when  compared  to  the 
effect  of  the  laser  gyro  errors,  the  results  of  the  sensi¬ 
tivity  analysis  on  the  laser  gyro  errors  is  more  importance 
than  that  of  the  accelerometer.  The  sensitivity  analysis 
indicates  that  the  bias,  misalignment,  and  scale  factor  errors 
are  responsible  for  almost  all  of  the  laser  gyro  error  induced 
INS  errors.  Using  once  again  the  -  1  nm/hr  position  error 
and  ^  3  ft/sec  velocity  error  specification  it  is  found  that 
a  50#  reduction  in  the  laser  gyro  errors  would  resutl  in  an 
INS  performance  improvement  of  50f°  which  would  be  enough  to 
meet  the  specification  above. 

A  final  overall  conclusion  is  that  although  the  high 
dynamic  environment  of  a  fighter-type  aircraft  causes  unac¬ 
ceptable  degradations  in  the  INS  performace,  it  does  seem 
feasible  that  there  is  a  bright  future  for  the  strapped-down 
INS. 
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In  performing  the  research  for  this  paper,  several 
recommendations  for  future  study  have  come  to  mind.  Among 
these  are  to  investigate  both  sensor  location  and  configu¬ 
ration.  The  sensor  location  is  important  because  of  the 
detrimental  effect  that  the  structural  modes  have  on  the 
INS  performance.  Also,  a  different  sensor  configuration 
could  reduce  the  magnitude  of  the  misalignment  error  effects 
on  the  INS  errors. 

Another  possibility  for  future  study  would  be  in  re¬ 
gards  to  the  simulation  used  in  the  present  research.  The 
present  simulation  model  is  an  F-4  aircraft.  It  would  be 
better,  though,  if  the  simulation  could  be  updated  such 
that  it  incorporates  and  F-15  or  F-l6  aircraft  model. 

A  final  recommendation  should  be  to  continue  to  refine 
and  improve  upon  inertial  sensor  technology.  The  sensiti¬ 
vity  analysis  indicated  that  a  big«step  forward  in  improving 
the  strapped-down  performance.  In  this  area,  however,  the  em¬ 
phasis  should  be  on  the  laser  gyro  rather  than  the  accelero¬ 
meter.  Also,  a  more  detailed  sensitivity  analysis  may  lend 
more  insight  into  the  minimum  reduction  in  sensor  errors 
that  is  necessary  to  meet  the  specification  used  in  this 
study. 
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Appendix  A 


Notation 

A  brief  description  of  the  notation  used  in  this  paper 
is  presented  here  as  an  aid  to  someone  unfamiliar  with  this 
notation.  Several  different  coordinate  frames  have  been 
used  in  this  paper.  These  coordinate  frames  are  indicated 
by  a  lower  case  letter  used  as  a  subscript  or  a  superscript 
with  the  superscript  indicating  the  frame  in  which  the  quan 
tity  is  coordinitized  while  the  subscript  indicates  the 
frame  the  quantity  is  referenced  to.  The  coordinate  frames 
used  are  as  follows : 

Inertial  frame  (i ) *  This  frame  has 
origin  at  the  earth's  center  and  is 
nonrotating  with  respect  to  inertial 
space. 

Earth  frame  (e)j  This  frame  is  non¬ 
rotating  with  respect  to  the  earth  and 
has  its  origin  at  the  center  of  the 
earth.  The  unit  vectors  are  defined 
as  follows:  xQ  points  toward  the  e- 
quator  and  the  Geenwich  meridian, 
y  pounts  toward  the  equator  and  90° 
e,  ze  points  toward  the  North  Pole. 

Computational  frame  (c)i  This  is  the 
frame  in  which  the  navigation  calcula¬ 
tions  are  done.  It's  origin  is  at  the 
center  of  gravity  of  the  aircraft  and 
the  unit  vectors  are  defined  as  fol¬ 
lows:  x  is  positive  up  with  respect 

V 


to  Earth's  surface,  y_  is  positive 
East  with  respect  to  Earth's  surface, 
zc  is  positive  at  some  angle  - 
(wander  angle),  from  Worth  with  res¬ 
pect  to  the  Earth's  surface. 

Body  frame  (b)i  The  body  frame  has 
its  origin  at  the  center  of  gravity 
of  the  aircraft  axes.  Its  unit  vec¬ 
tor  are  defined  along  the  aircraft 
axes.  Its  unit  vectors  are  defined 
along  the  aircraft  axes,  x^  is  out 
of  the  top,  y^  is  out  the  right  side, 
and  zfc  is  out  the  nose  of  the  aircraft. 

Vectors  are  denoted  by  a  line  underneath  the  vector 
symbol;  and  if  needed  have  a  superscript  to  indicate  the 
coordinate  frame  in  which  the  vector  is  expressed.  As 
an  example,  fc  is  the  specific  force  vector  coordinatized 
in  the  computational  frame.  When  angular  velocities  are 
used,  a  double  subscript  is  used  to  indicate  the  two  refer¬ 
ence  frames  between  which  the  angular  velocity  exists. 

&Le  is  the  angular  velocity  vector  of  frame  (e)  with 
respect  to  frame  (i). 

When  the  individual  components  of  a  vector  are  used, 
they  will  be  indicated  by  the  same  symbol  as  the  vector 
quantity,  but  without  the  underscoring  and  with  an  addi¬ 
tional  subscript  x,  y,  or  z  to  indicate  which  component, 
therefore,  w?ex  is  the  x  component  of  angular  velocity  of 
the  earth  frame  with  respect  to  the  inertial  frame  as  co¬ 
ordinatized  in  the  computational  frame.  When  all  compo¬ 
nents  of  a  given  vector  are  being  defined,  the  following 
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configuration  will  be  useds 


(A-l) 


Direction  cosine  matrices  are  used  to  transform  a  vec¬ 
tor  from  one  reference  frame  to  another.  They  will  be  de¬ 
noted  by  a  capital  C,  with  a  superscript  denoting  the  frame 
the  vector  is  being  transformed  to  and  a  subscript  denoting 
the  reference  frame  the  vector  is  transformed  from.  Thus, 


c£ib 


(A-2) 


The  values  of  the  individual  elements  of  a  direction 

• :  ■  J  •  . 

cosine  matrix  are  determined  by  knowledge  of  their  ini¬ 
tial  condition  and  the  angular  rates  between  the  two  frames 
of  interest.  Thus,  if  w^  (coordinatised  in  either  frame) 
and  the  initial  conditions  of  are  known,  then  C^(t)  can 
be  completely  specified. 

Another  notation  that  has  been  employed  is  that  for  a 
skew  symmetric  matrix.  When  a  cross  product  between  two 
vectors  is  desired,  x  n)e,  and  the  cross  product  is 
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coordinatized  in  a  single  frame ,  it  is  defined  as. 


rfek  e 
*ie  E 


(wie  *  £) 


e 


Therefore  W?*  is  defined  as, 
16 
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ek 

ie 
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iez 
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iez 
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we 

wiey 


-w 


rex 


-w 
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ley  lex 


(A-4) 


(A-5) 


A  final  note  on  notation  is  in  regards  to  the  time 
derivatives  of  vector  quantities.  Two  notations  will  be 
used.  These  are  dot  above  the  vector  or  a  lower  case  p 
preceding  the  vector. 


£c  *  fc  =  p£c  (A-6) 

When  the  operator  p  is  used,  if  the  vector  quantity  has 
a  superscript  to  indicate  the  frame  it  is  coordinatized 
in,  p  will  indicate  the  time  derivative  of  that  quantity 
with  respect  to  the  reference  frame  indicated  by  the  su¬ 
perscript.  If  no  superscript  is  used,  the  frame  that  the 
time  derivative  is  taken  in  will  be  indicated  by  a  sub¬ 
script  on  the  p.  Therefore,  the  Theorem  of  Coriolis  can 
t>e  stated  for  any  vector  V  as, 


P iv  B  Pb£  +  aib  x  V 


(A-7) 
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Derivative  of  the  Velocity  Differential  Equation 

The  velocity  differential  equation  used  in  the  strapped- 
down  local  vertical  wander  azimuth  navigation  model  is  de¬ 
rived  here  for  completeness.  This  equation  is, 

V  =  £  ♦  £  -  (2w?e  «■  wec)  x  V  (B-l ) 

From  Newtonian  mechanics,  the  second  derivative  of 
position  as  seen  from  an  inertial  frame  is, 

Pi  £  *  G  ♦  1  (B-2 ) 

where  f  is  the  specific  force  and  G  is  the  mass  attraction 
vector.  Defining  g  as  the  gravity  vector  composed  of  mass 
attraction  and  centripetal  components, 

£  =  £  -  wie  x  (wie  x  £>  (B-3) 

equation  (B-2)  can  be  written  as, 

P A2  E  =  £  ♦  1  *  SLie  x  (wie  x  £)  (B-4) 


Now,  defining  velocity  as  the  time  rate  of  change  of  the 
position  vector  as  seen  from  the  earth  frame. 


(B-5) 


Appendix  C 

Coordinate  Frame  Relationships 

Two  transformation  matrices  are  used  in  this  report. 
These  matrices  are  c£,  and  the  body-to-coraputational  trans¬ 
formation  matrix,  and  C^,  the  earth-to-computational  frame 

v 

transformation  matrix.  This  section  provides  the  deriva¬ 
tion  of  these  matrices.  The  earth-to-computational  frame 
transformation  is  considered  first. 

Figure  C-l  shows  a  diagram  of  the  relationship  between 
the  earth  frame  and  the  computational  frame.  These  two 
frames  are  aligned  when  ,  l,  and  1  are  all  zero.  At  this 
orientation  the  respective  x,  y,  and  z  components  of  each 
frame  are  aligned.  To  get  the  transformation  matrix  between 
these  two  frames,  the  directional  cosine  matrix  for  a  rota¬ 
tion  of  about  x  ,  of  about  y  ,  and  of  about  z  ,  (  u  =  wan- 

v  v  V, 

der  angle,  L=latitude,  l=longitude ) ,  is  found.  Now  defin¬ 
ing  the  order  of  rotations  as  first  the  z  axis  rotation, 

c 

then  the  y  axis  rotation,  and  finally  the  rotation  about 

the  x„  axis,  the  direction  cosine  matrices  are  multiplied 
c 

together  to  form  C^. 


10  0 

cosL  0  sinL 

cosl 

sinl  0 

Ce- 

0  coscx  sin 

0  -sinot,  cos 

0  1  0 

-sinL  0  cosL 

-sinl 

0 

cosl  0 

0  1 
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Defining  the  elements  of  c£ 

as 

CCE11 

CCE’12. 

CCE13~ 

Ce  * 

CCE21 

CCE22 

CCE23 

CCE31 

CCE32 

CCE33 

(C-2) 


results  in, 


CCEji=cosL  cosl 
CCEj2=cosij  sinl 
CCE1^=sinL 

CCE21  =  -sinc<  sinL  cosl  -  sinl  cose* 
CCE22=“Sin«  sinL  sinl  oosl  -cose* 
CCEg^sinoc  cosL 

CCE^1=-cosn  sinL  oosl  +  sintx  sinl 
CCE^2c-cos«<  sinL  sinl  -  sinoc  cosl 
CCE^^=cosck  cosL 


( C -3 ) 


From  this,  latitude,  longitude,  and  wunder  angle  can  be 
calculated. 


L  =  tan"1(CCE13  /(CCE^  +  CCE^)*) 
1  =  tan_1(CCEl2  /  CCE11) 

23  /  CEE33‘ 


tan_1(CCE5_  /  CCE__) 


(C-4) 

(C-5) 

(C-6) 


The  other  transformation  matrix  used  is  the  one  to 
transform  a  vector  from  the  body-to-the-computational 
frame.  Figure  C-2  shown  the  relationship  between  these 
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Figure  C-2, 


Wander  Azimuth  (computation)  and 
Body  Coordinate  Frames  (Kef  2*13) 
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frames.  As  was  done  with  the  earth-to-computational  frame 
transformation,  the  body-to-computational  frame  transforma¬ 


tional  matrix  is  defined  as  the  product  of  three  direction 
cosine  matrices.  Using  the  Euler  angles  for  the  respective 
rotations  results  in  a  rotation  of  (J)  about  z  ,  6  about  y  , 
and  finally  (  y  ♦  at  )  about  x  » 

v 


1  0 

0 

COS# 

0 

sin0 

cos  <P  -sind> 

0 

= 

0  cos(v"*-«) 

sin(v  +  <* ) 

0 

1 

0 

sin0  costf> 

0 

0  -sin(^+«c) 

cos (v'+a ) 

-sinP 

0 

COS0 

0  0 

1 

— 

- 

— 

- 

- 

- 

where  is 

the  roll  angle, 

is 

the  pitch  angle,  and 

(C-7) 


is  the  yaw  angle. 

Defining  the  elements  of  as, 


K 

CCB12 

CCli13 

=;ccb21 

ccb22 

ccb23 

£cb31 

ccb32 

CCB33 

results  in, 


(C-8) 


CCB^^  cose  cos  $ 

(C-9) 

CCB^2=-cose  sin4> 

(C-10) 

CCB^3=  sin  e 

(C-ll) 

CCB21=  cos(V+c0sin$>  -  sin(V+*Osine  cos$> 

(C-12) 

CCB22=  cos(y+rf)cos4>  ♦  sin(V+<A)sin  e  sin  4> 

(C-13) 

CCB23=  sin(^+^)cos  e 

(C-14) 

CCB3i=-sin(^+^)sin<t>  -  cos (V*/)sine  cos<|> 

(C-15) 

CCB32=-sin(V+^)cos $  ♦  cos(y^)sine  sin$ 

(C-16) 

CCB33=  cos(y+^)cos  e 

(C-17) 

153 


LOiL 


From  these  elements,  the  Euler  angles,  roll,  pitch,  and 
yaw  are  determined! 

$=  tan-1(-CCBl2  /  CCBn)  (C-l8) 
*  =  tan_1(CCB13  /(CCB^  ♦  CCB*2)*)  (C-19) 
V-  tan_1(CCB23  /  CCB33)  -  (C-20) 
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£0.  ^errors  in  position  and  velocity  to  each  inertial  sensor  is 
determined.  Using  this  sensitivity  analysis,  the  sensor 
parameters  required  to  obtairya'  performance  specification  of 
(3)1  nm/hr  position  error  and  %/J  ft/sec  velocity  error  are 
defined.  Finally,  recommendations  of  areas  for  improvement 
are  presented. 
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